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We study the singlet-triplet relaxation due to the spin-orbit coupling assisted by the electron- 
phonon scattering in two-electron SiGe/Si/SiGe double quantum dots in the presence of an external 
magnetic field in either Faraday or Voigt configuration. By explicitly including the electron-electron 
Coulomb interaction and the valley splitting induced by the interface scattering, we employ the 
exact-diagonalization method to obtain the energy spectra and the eigenstates. Then we calculate 
the relaxation rates with the Fermi golden rule. We find that the transition rates can be effectively 
tuned by varying the external magnetic field and the interdot distance. Especially in the vicinity of 
the anticrossing point, the transition rates show intriguing features. We also investigate the electric- 
field dependence of the transition rates, and find that the transition rates are almost independent 
of the electric field. This is of great importance in the spin manipulation since the lifetime remains 
almost the same during the change of the qubit configuration from (f, f) to (2,0) by the electric 
field. 

PACS numbers: 73.21. La, 71.70.Ej, 72.10.Di, 61.72.uf 



I. INTRODUCTION 



Spin-based qubits utilizing semiconductor quantum 
dots (QDs) are believed to be the prospective candi- 
date for quantum information processing^— Recently, 
silicon QDs have attracted much attention due to 
their outstanding spin-related properties^— Specifi- 
cally, the hyperfine interaction can be reduced by iso- 
topic purification^ The Dresselhaus spin-orbit coupling 
(SOC)2£ is absent thanks to the bulk -inversion symmetry 
and the SOC induced by the interface-inversion asym- 
metry (IIA) is rather weaker— Moreover, the electron- 
phonon interaction, which plays an important role in spin 
relaxation, is much weaker than that in III-V semicon- 
ductor QDs since there is no piezoelectric interaction in 
silicon. 14 All these special properties together lead to a 
long decoherence time in silicon QDs, which is of great 
help in the process of coherent manipulation and informa- 
tion storage. Furthermore, the physics in silicon is con- 
siderably rich owing to the presence of the valley degrees 
of freedom. Silicon has sixfold degenerate conduction 
band minima, which can be splitted by either strain or 
confinement in quantum wells into two parts: a fourfold- 
degenerate subspace of higher energy and a twofold one 
of lower energy. The twofold degeneracy can be further 
lifted by a valley-splitting energy due to the interface 
scattering. The valley-splitting energy has a strong de- 
pendence on the confinement length of the structurei 30 i 31 

Nowadays, spin qubits in silicon single and double 
QDs have been actively investigated^ - — In silicon single 
QDs, we have studied the singlet-triplet (ST) relaxation 
by explicitly including the electron-electron Coulomb 
interaction and the multivalley effect. Our results in 
the Voigt configuration agree quite well with the re- 
cent experiment by Xiao et al.. le Silicon double QDs, 
which have been proven very useful in exploiting the 



spin Coulomb blockade,^ have also attracted much at- 
tention. Recently, Raith et al£^ studied the magnetic- 
field and interdot-distance dependences of spin relaxation 
in single-electron Si/SiGe double QDs. Li et alM- cal- 
culated the exchange coupling between the unpolarized 
triplet and the singlet on the basis of a large valley split- 
ting. Culcer et alJ^- investigated the multivalley effect 
on the feasibility of initialization and manipulation of ST 
qubits, showing that the valley degree of freedom makes 
the physics of Si QDs quite different from that of single- 
valley ones. In their work, they analyzed the spectrum 
with the lowest few basis functions. As will be shown in 
this paper, these lowest basis functions are enough for the 
convergence of the energy spectrum under investigation, 
but are inadequate to study the ST relaxation time, simi- 
lar to the situation of III-V semiconductor-based QDs. 33 
The relaxation rates calculated with these lowest basis 
functions and the convergent ones differ by about four 
orders of magnitude. Therefore, it is necessary to employ 
the exact-diagonalization approach with a large number 
of basis functions in order to have the correct ST re- 
laxation rates. Moreover, the electron-electron Coulomb 
interaction, which is crucial to the energy spectra and 
the wavefunctions of the singlet and triplet eigenstates, 
was not explicitly calculated in the literature but rather 
given as a Hubbard parameter! 10 i 14 ' 23 

In this work, we calculate the two-electron ST relax- 
ation in SiGe/Si/SiGe double QDs by explicitly includ- 
ing the Coulomb interaction, the valley degree of free- 
dom as well as the source of the ST relaxation, i.e., the 
SOCs j 29 ' 34 We employ the exact-diagonalization method 
to obtain the energy spectrum and the Fermi golden rule 
to calculate the spin relaxation ratesi 33 i 35 Without los- 
ing generality, we focus on a large valley splitting case 
where the lowest singlet and three triplet states are all 
constructed by the lowest valley eigenstate. We investi- 
gate the double QD system with either a perpendicular 
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magnetic field (the Faraday configuration) or a parallel 
one (the Voigt configuration). We find that the energy 
levels of the lowest singlet and three triplet states have a 
strong dependence on the external magnetic field and the 
interdot distance. The perpendicular magnetic field af- 
fects the energy levels mainly by the orbital effect and the 
Zeeman splitting while the parallel magnetic field only in- 
fluence the energy levels via the Zeeman splitting due to 
the strong confinement along the growth direction. The 
interdot distance has a strong influence on the Coulomb 
interaction and the orbital energy. Besides, we also find 
that the transition rates of the channels among these four 
levels can be markedly modulated by the external mag- 
netic field and the interdot distance. Moreover, from 
the dependences of the energy spectrum on the magnetic 
field and the interdot distance, we observe the anticross- 
ing points between the singlet and one of the triplets. In 
the vicinity of the anticrossing point, the transition rates 
of the channels relevant to these two states present either 
a peak or a valley. Furthermore, we also study the effect 
of the electric field on the energy levels and the transition 
rates. With the increase of the electric field, we find that 
the configurations of the lowest four levels change from 
(1, 1) to (2, 0), where (n, m) indicates the numbers of oc- 
cupancy of the left and right dots. Very different from 
the magnetic-field and interdot-distance dependences, we 
find that the transition rates are almost independent of 
the electric field. This property is of great importance in 
the spin manipulation, since the lifetime remains almost 
unchanged during the variation of qubit configurations. 

This paper is organized as follows. We set up the model 
and lay out the formalism in Sec. II. In Sec. Ill, we em- 
ploy the exact-diagonalization method to calculate the 
energy spectrum and the Fermi golden rule to obtain the 
ST relaxation rates. We investigate the magnetic-field (in 
both the Faraday and Voigt configurations) , the interdot- 
distance and the electric-field dependences of the transi- 
tion rates. The features of the transition rates in the 
vicinity of the anticrossing points are also emphasized. 
Finally, we summarize in Sec. IV. 
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with m z representing the effective mass along the z di- 
rection. V(r) = V c + V z and P = p + (e/e)A = 
-ihV + (e/c)A with A = (- y B_ L ,xB_ L ,2yB ll )/2. H so 
stands for the SOC Hamiltonian, including both the 
Rashba 3 ^ and IIA2£~— terms. Then, one obtains 



H so = a Q (P x a y - P y <7 x ) + b Q (-P x a x 



PyOy) 
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where ao and bo represent the strengths of the Rashba 
and IIA terms, respectively. The Zeeman splitting is 
given by Hz = \g^B{B±a z + B\\a x ) with g being the 
Lande factor. H^ — eEx is the electric field term with 
an electric field applied along the x direction. Consider- 
ing that four in-plane valleys have much higher energies, 
we only need to include two out-of-plane valleys in the 
calculation. These two valleys lie at ±(/co) along the z 
axis with (k ) = 0.85(27r/asi). Here, as; = 5.43 A is the 
lattice constant of silicon. 10 H v in Eq. (p} describes the 
couplin g 30 ' 31 between these two valleys. For convenience, 
one uses the subscripts "z" and u z" to denote the valley 
at (feo) and the one at — (fco), respectively. 




FIG. 1: (Color online) Schematic of the double QDs. The 
black solid curve stands for the case without applied electric 
field. The red dotted curve represents the case with an electric 
field along the x direction. 2xq is the interdot distance. 



II. MODEL AND FORMALISM 



In our model, we choose the lateral confinement poten- 
tial as V c (x, y) = ±m t c^{min[(x - x ) 2 , (x + x ) 2 ] + y 2 } 
with rrit and coq representing the in-plane effective mass 
and the confining potential frequenc y 36 ' 37 A schematic 
of the double QDs is shown in Fig. [TJ The two dots 
are located at Rr : l = (±xo,0, 0) with 2xq being the 
interdot distance. Here R and L denote right and left, 
respectively. Along the growth direction [001], V z (z) is 
applied within the infinite-depth well potential approxi- 
mation. The single-electron Hamiltonian with magnetic 



To obtain the single-electron basis functions, we de- 



fine H = P * 2 t^ v — 1~ 5m — ^ ^( r )' The- ^-component of 
the Hamiltonian Hq can be solved analytically with the 
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eigenvalues being E nz = "g m x a 2 ■ Here a represents the 
half-well width. The corresponding eigenfunctions are 



frin^C* + a)], |z|< a 

otherwise 



(3) 



in which the index n z stands for the subband along the 
growth direction. In our calculation, only the first sub- 
band is included since the others have much higher ener- 
gies. It is very difficult to solve the in-plane part of Hq 
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analytically since the lateral confinement potential of the 
double QD, V c (x, y), lacks the symmetry of rotation. As 
the single-dot potential can be solved analytically , 17 ' 33 we 
solve the Hamiltonian of each dot separately instead to 
obtain the in-plane part of single-electron basis functions. 

We define Hl, 



Px +Py 



+ \m t ul{x±x ) 2 + \mtuly 1 
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The effective diameter can then be expressed as do = 
^hir/(m t uio). In the left dot, H L = P=c 2 ^ y — \-^m t ui 2 (x+ 
x ) 2 + \m t u 2 y 2 = ^-(p x , 2 +p' y 2 ) + \m t n 2 {x' 2 + y 2 ) + 
w B {x'p' y - yp x >), where x' — x + x , p' y = p y - w B x a m t , 
il = \Ju)q 2 + ojb 2 an d w s = eB±/(2m t ). We define 



h(P*' 2 + Py ) + §™t^V + V 2 )- One finds 



that the Hamiltonian Hl can be solved analytically in 
the system of polar coordinates while H' L can be solved 
analytically in the rectangular coordinate systemj 36 i 37 As 
it is much easier to calculate the Coulomb interaction nu- 
merically in the rectangular coordinate system and the 
term wb{x'p' — yp x ') can be treated perturbatively, we 
solve the Schrodinger equation of H' L instead of Hl to 
obtain the single-electron basis functions. One obtains 
the eigenvalue o 36 ' 37 

£ ni% = Ml(n x + n y + 1), (4) 

where n x ,y = 0, 1, 2, ... are the orbital quantum numbers 
of the x and y direction, respectively. The eigenfunctions 
are described as 
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W)]V /2 and 



ff nii „ s are the Hermite polynomi- 
als. Thus, the eigenfunctions in different valleys 
can be expressed as <f>%*fi nx = F^ ny (x,y)^(r) = 
F^ ny {.x,y)^ nz {z)e ±lk ° z u z ^(v) with u Z: - z (v) repre- 
senting the lattice-periodic Bloch functions^ Then 
one obtains a set of single-electron basis functions 
{{ffi!n v n J, wh ere are the eigen- 

functions in different valleys in the right dot and can be 
obtained by replacing L and xo in the left dot by R and 
— Xq. The orbital effect of the parallel magnetic field is 
negligible due to a strong confinement along the growth 
direction. 

In the present work, only H v is considered to con- 
tribute to the intervalley coupling since the overlap be- 
tween the wavefunctions in different valleys is negligi- 
bly small. 10 However, there still remain some controver- 
sies over the valley coupling nowadaysi 31 ' 38 ' 39 Here, we 
take <**f|#v|^> = and (*tf\H v \*tf) = 
according to Ref. |3ll Including this intervalley cou- 
pling, the single-electron eigenstates in the left dot be- 
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with V v representing the ratio of the valley coupling 
strength to the depth of quantum welli^I For the case 
of the right dot, one can get the corresponding single- 
electron eigenvalues and eigenfunctions by replacing L 
and xo in the left dot by R and ~xq. Then one 
obtains a new set of single-electron basis functions 



±L 



±R 



}}. It is noted that these basis func- 



tions are nonorthogonal and over-complete. 

Then we turn to the system of two-electron double 
QDs, where the total Hamiltonian is given by 



-fftot — {H e + H 2 + Hq) + Hp 



Hep 



H 



cp- 



(8) 



Here, the two electrons are denoted by "1" and "2". 

The electron-electron Coulomb interaction is given by 

2 

He = 47r£oK j ; ri _ r2 | with k standing for the relative static 

dielectric constant. H p = J2q\ f UjJ q\ a q\ a q\ represents 
the phonon Hamiltonian with A and q denoting the 
phonon mode and the momentum, respectively. The 
electron-phonon interaction Hamiltonian is described by 
#cp = EqA M qA«A + a- q A)e* qr and HI (i = 1,2) is 
given by Eq. (Q]). 

On the basis of the set of single-electron basis functions 
{{^nynJA^nynJh we construct the two-electron 
basis functions in the form of either singlet or triplet. For 
example, we use two single-electron spatial wavefunctions 
\n x m y in z in v ipi) and \n x2 n y2 n z2 n v2 p 2 ) (denoted as \Ni) 
and \N 2 ) for short; n v = ±; p = L/R) to obtain the sin- 
glet wavefunctions 



\s^) = (\n)-\m4P NlN2), 

' ' Kl ' ' " \\{\N l N 2 ) + \N 2 N 1 )), N X ^N 2 , 
and the triplet wavefunctions for N\ 7^ iV 2 



Nt =N 2 
^N 2 
(9) 
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= ^(|iViiV 2 ) - \N 2 N{)) ® (| tl) + I ItMll) 



|t! h) ) = 



-j={\N 1 N 2 )-\N 2 N 1 ))®\i\). (12) 



Here, the spatial wavefunctions of the first and second 
electrons in \NN') are denoted as N and N' in sequence. 
Specially, we denote p\ = p 2 = L/R as (2,0)/(0,2) and 
Pi Pi as (1, 1) configuration. The superscript (S) de- 
notes the valley configuration of each state. We define 
5 = ± for the valley indices of single electron states 
n„i — n v2 — ±, and 3 = m for n v i ^ n v2 . 

Then, one can calculate the matrix elements of two- 
electron Hamiltonian H^+H 2 + Hq in Eq. © and obtain 
the two-electron Hamiltonian matrix, where the matrix 
elements of the Coulomb interaction can be expressed by 
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(NiN 2 \H c \NM) = J2 ^ 1 ^^ l 4|0(<^,^ M> ^ B ^ l ^ B;]i ) 1 (13) 

7l,72,7j,72= 2: ! z 

where the superscripts 7$ and 7^ run over the two valleys, z and z, with jjj = 1 and jji = —77^ = 1. G is given by 

i- (6 luPl \e ik ' r \6 y ')' p ' 1 , , )(6 1 V P ' 2 , , le lk ' r l</> 72 ' P2 )* 

u ™ I l" s l^l'™.2n,2'i I 2'X 1 ii; i < 1 ''''n>^riy / £2 

(14) 



As two-electron basis functions are nonorthogonal, we 
also calculate the overlap between these basis functions. 
One finds that these two-electron basis functions can be 
divided into three independent subspaces according to 
the valley index, i.e., S = ± and m, as there is nearly 
no couping between them due to the negligibly small in- 
tervalley Coulomb interaction^ and overlap between the 
wave functions in different valleys J£ Then one can diag- 
onalize the eigen equation H^X — \S^ a 'X in each sub- 
space separately and obtain corresponding two-electron 
energy spectra and eigenfunctions, where and 
stand for the two-electron Hamiltonian and overlap ma- 
trix in the subspace with the valley index S (S = ± 
or m) respectively^ We identify a two-electron eigen- 
state as singlet (triplet) if its amplitude of the singlet 
(triplet) components is larger than 50 %. We use the 
similar way to identify a two-electron eigenstate as (2,0), 
(0,2) or (1,1) configuration according to the maximum 
amplitude. 

The transition rate from the state \i) to |/) due to 
the electron-phonon scattering is calculated by the Fermi 
golden rule, 

T ^f = YEl- M qA| 2 K/lxl<)r[fiqA«(6/-e i -na; qA ) 

qA 

+ (n qA + l)<5(e/-e l + fiw qA )], (15) 

in which x(q, ri, r 2 ) = e zqri + e zqr2 and n qA stands for 
the Bose distribution of phonons. In the calculation, the 
temperature is fixed at K and only the second term, i.e., 
the phonon-emission process, occurs. One finds that the 
transition between the eigenstates in different subspaces 
is almost forbidden because (/|x|i) m Eq. flB"]) is strongly 
suppressed due to a large intervalley wave vector (2ko), 
similar to the suppressions of the intervalley Coulomb 
interaction and overlap between the wave functions in 
different valleys mentioned above. 



III. NUMERICAL RESULTS 

From the Fermi golden rule [Eq. (|T5|) ]. one finds that 
the phonon energy is just the energy difference between 
the initial and final electron states. The energy difference 



studied here is much smaller than the energies of the in- 
tervalley acoustic phonon and the optical phonon4i Be- 
sides, the piezoelectric interaction is absent in silicon^ 4 
therefore one only needs to take into account the in- 
travalley electron-acoustic phonon scattering due to the 
deformation potential. In this work, both the TA and 
LA phonons are included. The corresponding matrix 
elements read M| intra Q = HD 2 Q 2 / (2d intra, q) with 
/3=LA/TA representing the LA/TA phonon mode. The 
deformation potentials for the LA and TA phonons are 
D LA = 3.93 eV and D TA = 2.48 eV, respectively^! 
The mass density of silicon d = 2.33 g/cm 3 42 The 
phonon energy f2^i n tra,Q = vpQ with sound velocities 
WLA = 9.01 x 10 5 cm/s and v T a = 5.23 x 10 5 cm/s4i The 
effective mass m t = 0.19mo and m z = 0.98mo with uiq 
being the free electron mass4£ The Lande factor g = 2, 44 
the ratio of the valley coupling strength to the depth of 
quantum well V v = 7.2 x 10 -11 m^ and the relative 
static dielectric constant k — 11. 9. 45 In the previous work 
on silicon double QDs^ i 14 ' 23 only the lowest few basis 
functions were included in the calculation. We find that 
these basis functions are enough for the convergence of 
the energy spectra, but inadequate in obtaining the cor- 
rect transition rates. The energy spectra calculated with 
the lowest few basis functions and the convergent ones 
differ by about 0.1 %, but the transition rates calculated 
with the lowest few basis functions differ by about four 
orders of magnitude from the convergent ones. Therefore, 
in our calculation, we employ the exact-diagonalization 
method with the lowest 1050 singlet and 3060 triplet ba- 
sis functions to ensure the convergence of the eigenstates 
and the transition rates. It is noted that all the ba- 
sis functions chosen are in the subspace with the valley 
index " — " since we focus on the case of a large val- 
ley splitting where a large effective diameter is taken to 
make sure that the lowest singlet and triplet states un- 
der investigation are in the subspace with the valley index 
"— " . This choice does not lead to the loss of generality as 
states with different valley indices are nearly decoupled 
as pointed out above. It is also noted that the Coulomb 
interaction was treated as a Hubbard parameter u in the 
previous works on silicon double QDs ) 10 ' 14 i 23 but in our 
work, it is explicitly calculated. One can obtain u from 
our calculation, e.g., u in the model by Culcer et al^ is 
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determined to be about 23 meV. 



A. PERPENDICULAR MAGNETIC-FIELD 
DEPENDENCE 

We first investigate the case of a perpendicular mag- 
netic field. We take 32 monoatomic layers of silicon along 
the growth direction of the quantum well, corresponding 
to the well width 2a = 4.344 nm. According to Eq. (0, 
a large valley splitting 2\Aj l \ = 0.83 meV is obtained. 
Then we choose the effective diameter do = 30 nm. With 
an electric field 30 kV/cm along the growth direction, 
one obtains the strength of the Rashba SOC induced by 
this electric field ao — —6.06 m/s and that of the IIA 
term b = -30.31 m/sM- 
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FIG. 2: (Color online) The lowest four energy levels vs. per- 
pendicular magnetic field B±_ in double QDs. The anticross- 
ing point between |S^ - ') and |T_ ') is shown and the range 
near this point is enlarged in the inset (the energies are sub- 
stracted by 48.89 meV). In the calculation, the interdot dis- 
tance 2xq — 20 nm. 

With the interdot distance 2xo = 20 nm, the lowest 
four levels are plotted in Fig. [2] as function of the per- 
pendicular magnetic field, denoted as IS^ - -*), |T+ ) (spin 

up), \Tq ) (zero spin) and \T_ ) (spin down) according 
to their major components. The energies of three triplet 
states are separated by the Zeeman splitting. From the 
figure, one finds that the singlet state intersects 



with the triplet states \T_ and |Tq ; ) in sequence 
with the increase of the magnetic field. The intersect- 
ing point between \S^) and |t1 _) ) (B± ~ 0.415 T) is 
an anticrossing point where there exists a small energy 
gap (~ 0.17 fieV) shown in the inset due to the Rashba 
SOC^i and the IIA term^ The intersecting point be- 
tween |£(->) and |T (_) ) (B ± ~ 0.98 T) is simply a cross- 
ing point. 

The anticrossing between IS^ - )) and |Tl ') can also be 
tuned by varying the interdot distance. We plot the en- 
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FIG. 3: (Color online) (a) The energy difference between the 
lowest singlet |5 ,( - _ '} and the lowest triplet \T_ ') vs. perpen- 
dicular magnetic field B± and half of the interdot distance 
xo in double QDs. The yellow solid cure is the position of 
the anticrossing between IS'' - '} and \T_ 

energy difference between IT"! - '} and \S^ 
difference of the Coulomb interaction u vs. half of the inter- 
dot distance xo at zero magnetic field and a magnetic field 
B± = 0.4 T. 



ergy difference between IS^ - ') and |Tl ^) as function of 
the magnetic field and the interdot distance in Fig.^a). 
In this figure, we also show the position of the anticross- 
ing between \T_ ) and |S ,< - - ' 1 ). It is seen that the mag- 
netic field where the anticrossing occurs decreases with 
the increase of the interdot distance. This can be un- 
derstood from the energy difference between |t£ ') and 
IS 1 ^'): A(Bj_,x ) = £ Q + u — gfisBj_, where £ is the 
orbital-energy difference and u comes from the contribu- 
tion of the Coulomb interaction. By solving the equa- 
tion A(B±,xo) = 0, one obtains the magnetic field B^ 
corresponding to the anticrossing point. To facilitate un- 
derstanding of the dependence of B^ on the interdot dis- 
tance, we plot the interdot-distance dependence of and 
u at zero and a specific magnetic fields in Fig.[3Jb). From 
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FIG. 4: (Color online) The lowest four energy levels vs. elec- 
tric field along the x direction. In the calculation, the interdot 
distance 2xo = 20 nm and the magnetic field B± — 1.4 T. 



trie field. Therefore, the electric field has the same effect 
on these states, leading to the energy differences among 
these levels insensitive to the electric field. 




B±(T) 



this figure, one finds that the orbital-energy difference 
decreases with increasing the interdot distance while u, 
which is insensitive to the magnetic field under investi- 
gation, shows an opposite behavior. It is noted that the 
increase of u is much smaller than the decrease of the 
orbital-energy difference £ Q . Therefore, with the increase 
of the interdot distance, the net contribution of £ Q + u 
decreases and correspondingly the magnetic field where 
the anticrossing occurs, i.e., B c ± , decreases too. 

In addition, the electric field can also effectively af- 
fect the energy levels. We apply an electric field along 
the x direction and plot the lowest four levels in Fig. @] 
One notices that the energy levels are weakly dependent 
on the electric field in the small electric field regime but 
show a rapid decrease when the electric field becomes 
strong. This electric field dependence agrees with that 
reported by Culcer et al^ qualitatively. Moreover, one 
also finds that the energy differences among these lev- 
els are almost independent of the electric field. These 
behaviors can be understood as follows. With the in- 
crease of the electric field, the configuration of these four 
states gradually changes from (1,1) to (2,0) according 
to their major components. In (1, 1) configuration, i.e., 
in the small electric field regime, the electric field sup- 
presses the single-electron energy in the left dot while 
raises it in the right dot. Therefore, the net contribution 
of the electric field is small and changes slowly with in- 
creasing electric field. However, when the electric field 
becomes strong, i.e., the states are in (2,0) configura- 
tion, the energy of each electron decreases while that of 
the two-electron Coulomb interaction increases with in- 
creasing electric field. It is noted that the increase of 
the energy of the Coulomb interaction is much smaller 
than the decrease of the energies induced by the electric 
field. As a result, the energy levels show a rapid decrease. 
Besides, we find that these four states always keep the 
same configurations regardless of the strength of the elec- 



FIG. 5: (Color online) Transition rates vs. perpendicular 
magnetic field. The inset zooms the range near the anticross- 
ing point of the channels S^ ) T\_~ ) and T < i~ ) T\r ) . In the 
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FIG. 6: (Color online) Transition rates vs. half of the interdot 
distance. In the calculation, the magnetic field B± — 0.3 T. 

We then calculate the ST relaxation rates together 
with the transition rates of the channels between two 
triplet states. From Fig. [SJ one finds that the transition 
rates can be markedly modulated by the magnetic field. 
In the vicinity of the anticrossing point (B± ~ 0.415 T), 
the transition rates show intriguing features. The rate 
between and \T_) shows a sharp decrease due to 

small phonon energy, which has been addressed in our 
previous investigations on GaAs and Si single QDsj 17 i 33 
The transition rates of other channels except the one 



between |7+ ') and |Tq ; ) present either a peak or a 
valley due to the large spin mixing between the singlet 



,(-) 
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1 5^ -*) and the triplet \T_ ), similar to the behavior we 
have investigated in single QDs. 17 Specifically, the transi- 
tion rate between IS^ - - 1 ) and (7+ ) and the one between 



|Tl } ) and |T(5 ') present a minimum while the transi- 
tion rate between IS^ - - 1 ) and |Tg - ^) and the one between 

\T_ ) and |Tl show a maximum. Far away from the 
anticrossing point, the variation of the transition rates 
can be well understood from the change of the phonon 
energyj 17 ' 33 i 46 

We also investigate the influence of the interdot dis- 
tance on the ST relaxation rates together with the tran- 
sition rates of the channels between two triplet states at 
B± = 0.3 T. The results are shown in Fig. [SJ We also 
find an anticrossing point between |5'^ _ - ) ) and at 
x ~ 12 nm. In the vicinity of this point, the behavior of 
the transition rates is similar to what we have obtained 
above by sweeping the magnetic field. The transition 
rate between \S^^) and \T_ ') is strongly suppressed and 
other transition rates relevant to these two states also 
show a rapid increase or decrease. Therefore, one can 
tune the ST relaxation by varying the interdot distance 
which can be controlled electrically in the experiment. 

Moreover, the electric-field dependence of the transi- 
tion rates is also studied. We find that the transition 
rates are almost independent of the electric field (not 
shown). This behavior can be understood from Fig. 0] 
where the energy differences among the lowest four lev- 
els are almost independent of the electric fiek h 17 ' 33 This 
property is of great importance in the spin manipula- 
tion, since the lifetime remains almost identical during 
the change of the qubit configuration from (1, 1) to (2, 0) 
by electric field. 
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FIG. 8: (Color online) (a) Transition rates vs. parallel mag- 
netic field By with the interdot distance 2xo = 20 nm. (b) 
Transition rates vs. half of the interdot distance with the 
magnetic field Bu = 0.3 T. 




FIG. 7: (Color online) The lowest four energy levels vs. par- 
allel magnetic field By in double QDs. In the calculation, the 
interdot distance 2xq = 20 nm. 



B. PARALLEL MAGNETIC-FIELD 
DEPENDENCE 



We also investigate the case with a parallel magnetic 
field along the x direction. The well width, effective di- 
ameter and strengths of both the Rashba SOC and the 
IIA term are chosen to be the same as the perpendic- 
ular magnetic-field case. We investigate the magnetic- 
field dependence of the energy spectrum in the absence 
of the applied electric field with the interdot distance 
2x = 20 nm. The lowest four levels are plotted in 
Fig. [71 We denote these states as l^")), |Tj -) ), |T ( ~ } ) 
and |T_ ) according to their major components. From 
the figure, one finds that \S^) and |T ( ] ) are almost in- 
dependent of the magnetic field while \t\t^} and |T_i ') 
show a linear dependence. This is because of the ab- 
sence of the orbital effect of the parallel magnetic field 
thanks to a strong confinement along the growth direc- 
tion. Therefore, the magnetic-field dependence is in- 
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volved only through the Zeeman splitting. For \S^ >) 

and |Tq the x-components of the total spin are al- 
most zero, leading to negligibly small Zeeman splitting. 
For |t| ') and \T^_ '), the x-components of the total 
spin are nearly ±1, which indicate that these two lev- 
els change linearly with the magnetic held. Besides, we 
also find an anticrossing point between IS*^) and |Ti ') 
at Sy - 0.54 T due to the SOCs. 

The influence of the magnetic held and interdot dis- 
tance on the ST relaxation rates is investigated. From 
Fig. [51 one finds that the behavior of the transition rates 
is quite similar to what obtained in the case of the per- 
pendicular magnetic field. Here, one also observes the 
anticrossing point between |S , ^~-') and |t1 ') by sweep- 
ing the parallel magnetic field and/or interdot distance. 
In the vicinity of the anticrossing point, the transition 
rates relevant to these two states also show a rapid in- 
crease or decrease. 



IV. SUMMARY 

In summary, we have investigated the ST relaxation in 
two-electron silicon double QDs with magnetic fields in 
both the Faraday and Voigt configurations. The electron- 
electron Coulomb interaction and the mutivalley effect 
are explicitly included. A large number of basis functions 
are utilized to converge the eigenstates and the transition 



rates. We find that the external magnetic field and the 
interdot distance have strong influence on the lowest four 
energy levels and consequently the transition rates can be 
effectively modulated by the external magnetic field and 
the interdot distance. Moreover, from the magnetic-field 
and interdot-distance dependences of the energy spec- 
trum, we observe ST anticrossing points. In the vicinity 
of the anticrossing point, a small energy gap exists be- 
tween the singlet and one of the triplet states due to the 
SOCs. The transition rates of the channels relevant to 
these two states show either a peak or a valley. Further- 
more, we also study the effect of the electric field on the 
energy spectra and the transition rates. We find that the 
configuration of the lowest four levels change from (1,1) 
to (2, 0) with the increase of the electric field. Differ- 
ing from the magnetic-field and interdot-distance depen- 
dences, the transition rates are nearly independent of the 
electric field. This is of great importance in the spin ma- 
nipulation since the lifetime remains almost unchanged 
during the manipulation of qubit configuration. 
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